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ABSTRACT 

Measurements of the ratio of the elastic electron- 

deuteron scattering cross section to the elastic electron- 

proton scattering cross section were made for low momentum 

transfers. Augmented by the work of Topping [5], these data 

2 -2 -2 

span a range of q from 0.05 fm to 0.50 fm . From these, 
the electric form factor of the neutron, G„ , was extracted 

hi 

n 

using three Lomon-Feshbach wave functions for the deuteron, 

2 

(LFl, LF5, and LF15) . Values of the slope of versus q 

■^n 

were measured to be 0.0218 ± 0.0022 for LFl, 0.0220 ± 0.0022 
for LF5, and 0.0185 ± 0.0022 for LF15. Although the data 
seems to prefer model LF15, all were found to be in agree- 
ment with the thermal neutron-electron slope of 0.0193 ± 
0.0004. The root mean square structure radius of the 
deuteron was also extracted, using a method devdsed by 
Schumacher, and was found to be 1.9665 ± 0.0045 fm. Rela- 
tivistic corrections proposed by Casper and Gross [6, 18] 
were compared with the relativistic correction of Friar [8] 
and no significant differences were found in the region 
investigated. 
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I. 



INTRODUCTION 



The deuteron being composed of one proton and one 

neutron is the only stable two body nucleus and, as such, is 

the simplest of the nuclei except for the hydrogen nucleus. 

The proton exists as a stable particle in the form of a 

hydrogen nucleus, but the isolated neutron is unstable. In 

order to obtain information concerning the neutron, this 

work will compare electron scattering from the deuteron and 

proton and extract the electric form factors of the neutron, 

Gg , which gives information on the charge structure of the 
n 

neutron. 

o 

This experiment measured at various values of q" 

n 

(the square of the Lorentz invariant four-momentum trans- 
ferred from the scattered electron to the target nucleus) 

-2 -2 

ranging from 0 . 4 f m to 0 . 5 fm , in order to resolve 

differences among several earlier experiments [1,2, 3, 4]. 

The results of this work augmented with those of Topping [5] 

provide values of the neutron charge form factor for the 

2 - 2-2 

momentum range q = 0.05 fm to 0.50 fm . By comparing 
these G_ , which are obtained with the use of several dif- 

Xj 

n 

ferent wave functions, to the G„ predicted by the very 

hi 

n 

accurate thermal neutron experiment of Krohn and Ringo [4], 
a model that best represents the structure of the deuteron 
may be found. Three different models containing D-state 
contributions of 4.57, 5.20 and 7.55% were tested and 
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although the results are not conclusive, the data favors 
the model with 7.55% D-state. Two different relativistic 
corrections were investigated, the first developed by 
Gross [6] and modified by Schumacher [7] and the second 
developed by Friar [8]. Results using both methods are 
reported here. For the momentum transfers of this experi- 
ment, it is shov/n that the differences between the two 
methods are negligible. From the slopes of G_ for the 

Ij 

n 

various deuteron models used, a mean structure radius of 
the deuteron is extracted employing a method developed by 
Schumacher [7] . 

The organization of this thesis is as follows; the 
concepts of electron scattering and the definitions of 
terms needed to describe the scattering from protons and 
deuterons are contained in Chapter II. 

Chapter III discusses briefly the experimental equipment 
and the procedures followed in order to obtain -the data. 

The contents of Chapter IV discusses the methods employed 
for data reduction. Included are discussions of corrections 
and normalizations made to the data, radiative and ioniza- 
tion corrections to the scattering spectra, and a detailed 
analysis of the sources of errors. 

Chapter V presents the results using several different 

wave functions and both methods for incorporating the 

2 

relativistic correction. Plots for G versus q are 

£ 

n 

included for all cases. Also included is a Rosenbluth Plot 
as an experimental verification of the conventional scaling 
law assumptions of Chapter II. 
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Chapter VI presents the conclusions of this work and 
recommendations by the authors for future experiments. Five 
appendices amplify details of several topics discussed in 
this thesis. The last appendix is a compilation of the 
FORTRAN codes used in the data analysis. 
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II. THEORY 



A. NEUTRON CHARGE FORM FACTOR 

The cross section for an extremely relativistic electron 
scattering from a spinless point nucleus with no magnetic 
moment has been calculated by Mott [9] f 



aM(e) = 



fda\ 


(ze^l 


1 dft .. 


2E . 


V /Mott 


1 



2 e 

cos 2 1 

.46 C 
sin 2 



(II-l) 



where 



2E. , 

C = 1 + — ^ sin^ 

M c 
n 



is the recoil factor. 

The experimental cross section is related to the Mott 
cross section by 



'd£\ ^ 



dn j 



exp 



d^J 



G^(q^) , 



Mott 



where G(q'^) is the total form factor of the nucleus, which 
in the case of charge scatterings is a function only of the 
square of the four-momentum transfer. The form factor takes 
into account the finite extent of the nucleus . 

A theoretical expression for the experimental cross 
section for a proton has been developed by Rosenbluth [10] 
and is given in terms of the Mott cross section and total 



form factor G^ as developed by Hand^ejb al . [11] 
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(II-2) 






where a^(0) is the proton Mott cross section and 
M 



2 _E 



G = 
P 



1 + T 



^ + 2t tan^ T '/ 

2 M 



where and Gj^ are the electric and magnetic form factors 

P P „2 

of the proton respectively and t = . 

4M 

P 

In principle. Equations II-l and II-2 may be used to 

measure G^ and G^ by doing experiments at constant q with 
P P 

various values of incident energy and scattering angle. 
However, for low values of q, the magnetic term is small (at 
most 15% of the total) so that it is more convenient to cal- 
culate the magnetic scattering using the scaling law developed 
by Lehmann et al . [12] and others [11] , which was shown to 

be valid at low momentum transfers by Brooks and Sessler [13] . 

In terms of the proton magnetic moment, y , the' scaling law 

P 

is given by 



3 ., — y G_ 

M ^p E 
P ^ P 



(II-3) 



Utilizing Equation II-3, Equation II-2 can be simplified 



to 



E 

Gp = [1 + Typ {1 + 2(1 +t) tan^ |} ] 



(II-4) 



. 2 

Since the xy^ term in the brackets is due only to magnetic 

2 2 6 

factors, defining = xyp{ 1+2 (1+x) tan , equation II-4 
can be written 
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S = 1+^ ' 



Rewriting this and defining a proton magnetic correction 
term, = y^I — , an expression for is found. 



gI = rP 

E p M 

P ^ 



(II-5) 



From Equation II-2, it can be seen that G^ can be 
obtained from the ratio of the experimental cross section 



to the Mott cross section and, G„ can be written in terms 

hi 

p 

of measurable quantities. 

daP 






dfi J.P 

p / ,^ X M 



(II-6) 



An analogous method is used to develop an expression for 

2 

G_ , the electric form factor of the deuteron. Jankus has 
d 

developed an expression for the cross section for electron 
scattering from the deuteron [14]. Written in a form that 
includes form factors to account for the finite extent of 
the nucleons, Jankus' results [15] are 



S- = <’^(9)^ (q^) +1 (q2) +1 (q 2 )[l+ 2 (l+n)tan 2 |]), 

a d d 

(II-7) 

2 

where r\ = — and (q^) ^ (q^) / Gj^ (<3^) respectively 



4M 



the charge, quadrupole and magnetic moment form factors of 
the deuteron. They are given by; 



= (G^ + G^ ) C 



E 



E 



n 



'E 



= (G^ + G„ ) C 



E 



n 
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r 







2 




M, , 

fr <=E + S ' =M > =S • 

P P n p n ^ 



Here, C„ is the charge structure factor, is the quadru- 
t U 

pole structure factor, accounts for the magnetic contri- 
bution to scattering arising from charge convection in the 
deuteron and Cg explains the contribution to scattering by 
the magnetic moments of the proton and neutron. These four 
factors C„, C_, C- , and C„ are all model dependent in that 

III U ±J O 

they arise from integrals involving S and D state wave 
functions of the deuteron [5] . Values for C„, C^, C,. , and 

hi U Li 

Cg for momentum transfers studied were supplied by Lomon [16] 
for the Lomon- Feshbach Models LFl, LF5 and LF15 [17] which 
contain percent D-state contributions of 4.57, 5.20, and 
7.55 respectively. 

If u(r) and w(r) are the normalized S and D-state wave 
functions, then Cg is given by 

OO 

= J [u^(r) + w^(r)] (3^) dr , 

o 

where is a spherical Bessel function, corresponding expres- 
sions for Cq, Cg, and Cg may be found in Ref. [39]. In 
the static limit (q=0) , the above functions are given by 

= 1 , Gq^(O) = 

Cg(0) = 1, Cq(0) = m2q , 

Cg(0) = 1 - 3/2 Pg , Cg(0) = 3/2 Pg , 
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where Q is the quadrupole moment of the deuteron, y and y 
are the proton and neutron magnetic moments respectively, 

oo 

2 

and = / w (r)dr is the D-state probability of the 
o 

deuteron . 

Relativistic corrections to the deuteron scattering 
problem have been calculated by Gross [6] , and Casper and 
Gross [18] and Friar [8]. These corrections account for 
relativistic modifications of the deuteron wave functions 
and the nucleon current incorporating not only Lorentz 
contraction in two dimensions, but also the fact that the 
equal-time wave function in the rest frame is not the same 
as the equal time wave function in the moving frame. 

Schumacher and Bethe [7] and Elias et al. [15] have shown 
that the relativistic corrections of Casper and Gross [6,18] 



to the deuteron form factors, to order — ^ , enter as a 

, 2 '’«P^ 

factor (l--j) on the right hand side of Equation II-7. 

Friar has shown that the relativistic correction enters 

in two ways; first, a factor (1 +t) ^ should multiply the 

right hand side of Equation II-7; second, C„ must be eval- 

rj 

2 * 2 

uated at an effective q = q /(l+n). A comparison of the 
results using both formulations may be found in Section V. 

In Appendix D, it is shown that by expanding Friar's 
correction in a binomial series for low q , Schumacher's 

4 

correction differs by less than 4 parts in 10 for the 

. 2 

highest q investigated. Consequently, the change in C^ is 
negligible at low q . The relativistic corrections of 
Casper and Gross [18] in the form developed by Schumacher 
[7] are used in this work. 
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When Equation II-7 is used to analyze the deuteron 

experiments, a scaling law for G is used [3]. 

^d 

Md 

^ ^d ^E , 

d P d 

The quadrupole term is dropped since it is less than 5 x 10 

2 

in comparison to the total at the highest value of q 
reached in this experiment [38] . Equation II-7 may be 
written in the form 



da 

dQ 



= I ‘■d GL^[l+2(l+n) tan^ |]j(l - |) . 



d d 

2 2 r 2 0 T 

Defining ^+2 ( 1+n ) tan and noting that 

d 



d^ - ^d ' 



where is the total form factor of the deuteron, then 



4 = * =E > ‘^E <1 + '^d> <1 - ¥ 

p n 



(II-8) 



Additionally, by defining as the deuteron magnetic and 

relativistic correction factor, the following is obtained: 

1 f 1 y 



K 



M 1+C, 

r d 



From Equation II-8, the deuteron electric form factor, G„ , 

d 

is obtained in terms of ; 

M 

r 



da 

^2 _ ^2 ^,d _ dQ 

E , d • M d , » . L-x 

d r . a., (6) r 



M 



‘'m = <°E * °E =E • <”-S> 

p n 
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The ratio of the proton and deuteron electric form fac- 



tors, is expressed in terms of G^, and as 

n p 

=E<^d> (=E <<3^ * =E <<3d> ' , 

d p n ^2 



=E «3;> 



(II-IO) 



P 



Combining Equations II-6, and II-9, with Equation II-IO, 
the ratio in Equation II-IO can be written in terms of 
measured quantities only. 



E n 

p ^ ■ 



^ 

do 


^5 


r daP n 


d d« 




T,P dfi 


M d,„. 




t- 

Q 

CD 

1 



(II-ll) 



It is important to note that q^ and q^, which are the 
four momenta transferred to the deuteron and proton respec- 
tively have different values. The incident energy was 
adjusted to make the final electron energies from both the 
deuteron and proton scatterings equal. This insures that 
the transmission efficiency of the spectrometer and the 
detector efficiencies were constant for both gases, but the 
four momentum transfers were slightly different. The 
deVries b' fit [19] for the proton electric form factor was 
used to adjust for this slight difference in momentum 
transfers . 



=E <<3d> 



=E ‘‘Jp' 
P ^ 





deVries 



(11-12) 



The ratio in Equation II-ll may be written 
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I 

i 

I 

I 

( 

i 



Ge <‘3d> + Gg '^d> 
_JB 2 

P 



do^ P/Q\ K.. 

dR- Jir 

'E I 0^(6) kP 



Ge %> 

P 

Ge •‘5d> 



deVries 



Solving this for (q^) yields: 

n 

i 

.d 



=E <Gd> = Gg <<5d> 
n p 



exp 

■^Mott rP 




Ge <'3^ 

p 

Ge '<3d> 

p 



1.0 



(11-13) 



deVries 



^exp experimental ratio of the deuteron cross 

section to the proton cross section, and R is the ratio 

of to 0^(0). Note that o^(0) and were calculated 

at the incident energy of the electron on the hydrogen, which 
was slightly higher than the energy of the electrons on the 
deuterium. 

Writing 




R 

exp 

^Mott 







9 

■ deVries 



Equation 11-13 can be rewritten in the more convenient form. 



Ge <Gd' 




1.0 





(11-14) 



2 

In this form G„ (q,) includes the effects of the Gross rela- 
tivistic corrections to the deuteron form factor. 



B. NEUTRON ROOT MEAN SQUARED CHARGE RADIUS 

In the Breit or "brickwall" frame, the form factors 

Gg and G^^ can be separated. The Breit frame is the 
n n 
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reference frame in which the incident and final electron 

energies are equal, and the initial and final momenta of 

the deuteron are equal but' oppositely directed [20]. 

2 

In this frame (q ) is the Fourier transform of the 

n 

charge distribution [21,22] and can be written as 



(q^) =/" p(r) . 

n o 



The exponential term may be expanded in a power series 

2 n th 

about q = 0 and the coefficients of q are the n moments 

of the charge distribution. For symmetric charge distribu- 
tions, only the even moments are non-zero. is given by 

n 



1 . 2 . . 2 . 



G (q ) = G (0) - ^ <r > q + 
n n 



120 



+ . . . , (11-15) 



where is the mean square charge radius of the charge 

. . . 2 

distribution. The slope of G„ as a function of q for low 

Xli 

n 

momentum transfers is related to the mean square neutron 
charge radius by 



G 

dq n 



6 ^CH 
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Ill . EXPERIMENTAL EQUIPMENT AND PROCEDURE 



A. NAVAL POSTGRADUATE SCHOOL ELECTRON LINEAR ACCELERATOR 

The. 100 MeV NPS LINAC is an electron accelerator similar 

in design to the Stanford University Mark III linear accel- 
erator described by Chodorow et al . [23] . The thirty foot 

accelerator is powered by three Klystrons, each of which 
receives up to twenty-two megawatts of peak power. The 
beam energy may be varied from 30 to 105 MeV with an average 
maximum current of about lyA. The machine is pulsed at a 
rate of 60 pulses per second, and each pulse is about one 
microsecond in duration. 

The beam energy is selected by the beam deflection 
magnets and energy defining slits. The field of these mag- 
nets is measured by a nuclear magnetic resonance (NMR) probe. 
With the aid of quadrupole magnets, the beam is- focused to 
a beam diameter of about 2 mm at the target. A detailed 
description of the NPS LINAC is given by Barnett and 
Cunneen [24] . 

B. SPECTROMETER AND COUNTING SYSTEM 

The spectrometer, described in detail by Oberdier [25] , 
is a movable 120®, 16 inch radius, double-focusing magnetic 
spectrometer able to detect electrons scattered at angles 
between 45® and 150® in steps of 15®. Ten NE102 plastic 
scintillation counters span part of the focal plane of the 
spectrometer, where each detector serves as the source of a 



20 




I 



counting channel. The width of the counters allows a momen- 
tum spread, Ap/p, of about 0.3% for each channel. Two large 
detectors, spanning the same field as the ten smaller 
detectors, served as dual coincidence backing counters. 

Thus an electron passing through both backing detectors and 
a front detector will register as one count on the Jorway 
Model 1836 Dual Channel Scaler located in the control room. 

A triple coincidence pulse will register on a scaler only 
if it occurs during the time that the counter gating circuit, 
gated with the accelerator gun pulse, is open. This pre- 
vents background counts from registering between electron- 
beam pulses. 

Scalers are capable of recording periodic pulses up to 
150 MHz, but due to the pulsed nature of the beam, the 
count rate was held below 60 counts/sec on all runs to avoid 
large count rate losses. The count rate correction was 
found to be linear up to 100 counts/sec and is applied in 
the codes BERHFORM and TRADFORM. It was decided to use only 
counters 4 through 7 in order to consistently follow 
Topping's [5] procedure since the data points accumulated by 
Topping were used in the present work. Sufficient data were 
taken such that the ratios found from channels 4 through 7 
were well defined. 

The gas target cell used was a two inch diameter steel 
cylinder with 1.0 mil Havar windows positioned to allow the 
incident electrons to enter and the scattered electrons to 
exit. Both the H 2 and D 2 gases were maintained at liquid 
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nitrogen temperature and about 150 psig. The gas target 
system is described in detail by Topping [5] . 

C. DATA ACCUMULATION 

Before entering into the actual procedure for obtaining 
data several terms must be defined. A "point" is the set 
of data obtained at one spectrometer setting for either 
the proton or deuteron. It consists of the number of 
counts registered in each of the ten counters, the integra- 
tion capacitance and voltage, machine and spectrometer 
energies and the time required to take the data. A "peak" 
is the collection of data associated with one of the ten 
channels, e.g., "peak 4" refers to the data collected by 
counter number 4. A "run" is the collection of all the 

points of all the peaks for both the deuterium and hydrogen 

2 

at one value of q and scattering angle. 

2 

Topping [5] reported results of 21 runs from q = 0.05 
-2 -2 

fm to 0.35 fm . The results of his runs and the 12 runs 

made by the authors will be combined to extend the data set 

-2 

to .05 fm . The data collected by the authors were mainly 
obtained under the same conditions and same methods as used 
by Topping [5]. 

All data for each point were recorded automatically (see 
Section IV-A) . A run consisted of between 50 and 75 points 
first taken with hydrogen in the target and then repeated 
with deuterium. 

The spectrometer was first "saturated" at a field value 
corresponding to an electron energy of 100 MeV before data 
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were taken on either or D 2 to insure that the same 
hysteresis loop was followed on both spectra. The spectro- 
meter was then set such that the elastic scattering spectrum 
was well above channel 4. The energy selected by the 
spectrometer was stepped down 0.05 MeV between each point. 
This process was continued until the energy selected by 
the spectrometer was about two full widths at half maximum 
below peak 7. Both the and D 2 spectra were obtained in 
the same manner. 

Two additional bits of data were recorded by hand, the 
Secondary Emission Monitor (SEM) efficiency and the pressure 
in the gas target. The SEM efficiency was measured when 
each of the peaks registered a maximum number of counts; 
therefore, a total of eight values were taken for each run 
and they varied less than 0.1%. 

The above data were always taken under conditions that 

were as similar as possible for the deuteron and the proton. 

2 -2 -2 

The authors made runs from q = 0 - 4 fm to 0 . 5 fm . The 
time required to complete a run in this experiment was 
greater than for the experiment of Topping because of the 
decrease of the cross section with q , and the fact that 
back angles had to be used exclusively. Consequently H 2 and 
D 2 spectra were accumulated on two successive days. This 
was a deviation from Topping's procedure, but this method 
gave data that reproduced within statistics. 

The beam charge integration was always set such that the 
hydrogen peaks contained from 10,000 to 20,000 counts in a 
single channel. 
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As in Ref . [ 5 ] , 

counter efficiencies 
for each channel. 



this work eliminates the need for 
by calculating a ratio of 



do's /doP 



/: 



dQ / dQ 
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IV. DATA ANALYSIS 



A. DATA REDUCTION 



The data collected on a Model 33 teletype was transcribed 
onto computer data cards compatible for use on the IBM 360-C. 
The data cards were organized such that there were two cards 
for every point obtained at a given spectrometer setting. 

The first card contained the counts for each of the ten coun- 
ters, and the second contained the number of background 
counts, the integration time, incident electron energy, the 
voltage integration, the energy setting for the spectrometer, 
the capacitor setting, and the integrated beam charge as 
measured by the SFM. The hydrogen and deuterium data were 
treated separately with FORTRAN programs BERHFORM and 
TRADFORM respectively which are listed in Appendix E. These 
programs are rivised versions of TOP#FORM [5]. They are used 
primarily to calculate the cross section measured by channels 
4 through 7, but they also make count rate corrections and 
charge monitor normalizations. Charge monitor normalization 
was required because the automatic termination of the inte- 
gration for each of the spectrometer settings did not always 
stop the counting system at the same value of beam charge. 

The correction was 



V 



N (normalized) = N (actual) 



reference 
V 



(IV-1) 



actual 

where is the actual voltage at the end of integration 

and is an arbitrary standard voltage used to nor- 

malize the hydrogen and deuterium counts. 
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These programs were also used to plot the normalized 
count rate corrected spectra on the CALCOMP plotters at the 
NPS Computer Center. For the first plotting run, certain 
additional data, which are necessary to calculate cross 
sections, (such as the radiation corrections and densities) 
were given arbitrary values. The count-rate corrected and 
normalized counts for both deuterium and hydrogen peaks 
were plotted. Using these plots, the position of each peak 
and the line shape of the electron spectrum was visually 
determined. A smooth curve was fit through the points and 
those points that were affected by energy instabilities were 
appropriately adjusted. Typical plots of both raw data and 
smoothed curves are reproduced in Figures IV-1, IV-2, IV-3 
and IV-4. 

The radiation correction program TOPRADCR [5] requires 
the following inputs: 1) the peak energies; 2) the pressure 

in the gas target for hydrogen and deuterium; and 3) AE which 
is the integration interval from the peak energy to the 
final spectrometer setting at which data was taken. These 
data were obtained from the before mentioned plots . TOPRADCR 
also used the composition and thickness of the gas target 
windows and the virial coefficients of hydrogen and deuterium 
as stored constants. TOPRADCR calculated the Schwinger 
radiation, bremss trahlung and Landau ionization corrections 
along with the density of the gas for each channel (peak) . 

The necessity for the radiation corrections will be discussed 
later in Section IV-C. Average background was determined 
from the spectrum above the elastic peak. 
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Figure IV-1 
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Figure lV-2 
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Figure IV- 3 
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The total radiation and ionization corrections, the 
gas density, the SEM efficiency and the background counts 
for each channel were typed on computer cards and were used 
as input data to BERUFORM and TRADFORM along with the 
appropriate data decks for hydrogen and deuterium. These 
programs then recalculated the elastic cross section and 
the statistical error in that cross section for each channel 
considered. The calculation of the cross section is 
discussed in Section IV-B, and the method for calculating 
statistical error is derived in Section IV-D. 

Through the use of BERHFORM and TRADFORM, four cross 

sections (one for each channel) were obtained for each of 

the hydrogen and deuterium peaks. These cross sections 

were then put in the form of a ratio (deuterium cross 

section divided by the hydrogen cross section), and a 

weighted mean ratio, R , was calculated. The electric form 

m 

factor of the neutron, G„ , for each run was de.termined 

n 

using R^^. The values of the proton electric form factors 
and the proton and deuteron Mott cross sections and magnetic 
corrections were obtained from programs TOPROTON and 
TOPDEUTR [5] . 



B. CALCULATION OF THE ELASTIC CROSS SECTION 

The experimental cross section for each channel is 



N K 
sc s 






K 



exp 



N. N^ Afi 
in t 



^ cm^ ster ^ , 



(IV-2) 
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where 



N is the number of scattered electrons counted 
in a particular channel, 

N. is the number of incident electrons, 
in 

2 

is the number of target atoms per cm , 

An is the solid angle defined by the spectrometer 
slits , 

K , K. , and K. are radiation and ionization cor- 
s' b 1 

rections (see Section IV-C) , 
is given by 

pt N 

N. = — , (IV-3) 

t A 

2 

where p is the density of the target (g/cm ) , 
t is the target thickness (cm) , 
is Avogadro's number, and 
A is the atomic weight of the target gas. 



N. is determined by 
in ^ 



N. 

in 



Q 

e 



CV 



e e 



(IV-4) 



sem 



where Q is the' total charge of incident electrons, 

-19 

e is the electron charge (1.60210 x 10 coulombs), 
V is the voltage of the integrator in volts , 

C is the capacitance of the integrator in farads , 

^sem^^ the SEM efficiency (about 5.8%). 

In this experiment a ratio of cross sections is measured, 
so that several of the parameters of Equations IV-2 and IV-3 
which are not absolutely determined, i.e., solid angle and 
target thickness, cancel with a subsequent elimination of 
systematic errors. 
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N is the total number of scattered electrons, 
sc 

Area under scattering peak 
_ (in counts - MeV) 



(IV-5) 



where E is the energy width of an individual detector. 

In order to obtain , the number of counts in a par- 
ticular channel must be integrated over the energy range 
of interest. 



n-1 

N = E 4 [ N.+N.^, ] (E. -E.^, ) 

sc 2 '• 1 1+1 1 1+1 



/o. 



0034877 E 



(IV-6) 



where is the number of counts in the channel for the i^^ 
spectrometer setting, and E^ is the energy of the i^^ spec- 
trometer setting. 

Combining Equations IV-3, IV-4 and IV-6 in Equation IV-2, 
the experimental cross section may be written as 



1 1 
2 ^^i‘*’^i+l^ ^^i“^i+l^ 0.00034 877E^(1.0+0. 0039 (n-5)) 



K K. K. 
s b 1 



exp 



CV 



ee 



sem 



Pt 

A 






(IV-7) 



where the expression 



E^ = E^ (1.0 + 0.0039 (n-5)) 



(IV-8) 



expresses the electron energy seen in counter n in terms of 
the spectrometer dial setting E^. Note that counter 5 sees 
the electron energy read from the spectrometer control dial. 

Count rate corrections were applied to N^ by using a 
linear correction empirically determined from a thick target 



33 



experiment. The empirical correction is 



, i (recorded) 

N (count rate corrected) = i.q - 0.00097 N . /t. ' 

i' 1 

where N^^/t^ is the average count rate. Background was deter- 
mined by averaging the counts for each peak above the start 
of the elastic peak. Charge integration normalization V7as 
done as given in Equation IV-1. 

is defined as the counts in the i^^ channel, corrected 
for background (B) , count rate and charge integration 
normalization . 



C. 

1 



V 

^i V 



ref 



act 



1 - 



0.00097 N./t 



- B 



(IV-9) 



The final expression for the cross section in terms of 
defined quantities is 



1 

— T 


(VCi^l) 


(E.-E.^,) K K,K. 
X x+1 s b X 




2 . , 
r. _ 1=1 


0.0034877 (E. /E . ) 

xn' X 




exp 


1 'I 




/ 




V ^semj 


A 





(IV-10) 



where Afl is the solid angle and has a value of 1.831 x 10 



steradians . 



C. RADIATION CORRECTION CALCULATIONS 

Both the proton and the deuteron cross section were 
derived with the assumption of one virtual photon exchange. 
This process is represented in Figure IV-5. For the deuteron 
the assumption of the impulse approximation is also made, as 
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represented in Figure IV-6. However, in the actual scatter- 
ing process, more processes than those illustrated in the 
two figures occur. The experimenter must correct for 
several competing processes in order to be able to compare 
experiment and theory. These competing processes are 

1. The emission of real photons and the reabsorption 
of virtual photons in the field of the scattering nucleus. 

The correction to the cross section for this process was 
first derived by Schwinger [26] and is therefore called the 
Schwinger Correction. 

2. Bremsstrahlung emission in the target before and 
after the nuclear scattering of interest. 

3. Excitation and ionization of material in the target. 

In order to obtain the cross section for a nuclear 

scattering event one would have to continue down to zero 
energy so as to count all the electrons which have some 
probability of losing all their energy in the three processes 
previously discussed. In practice the measurement of the 
cross section is stopped at some AE below the elastic peak, 
so it is necessary to include a correction to account for 
the undetected electrons with energy less than E^ - AE. 

Where E^ is the energy of the elastic peak and AE is the 
energy from the elastic peak to the energy corresponding to 
the last spectrometer setting, usually two to four full 
widths at half maximum. 
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C. RADIATION CORRECTION CALCULATIONS 



1 . Schwinger Correction 

The Schwinger correction was generally about 15% in 
this experiment and is the largest of the three corrections. 

The Schwinger correction accounts for processes 
which occur in addition to the one virtual photon exchange. 
Schwinger's original work was later improved by Tsai [27] 
and by Meister and Yennie [28] with the inclusion of the 
effects of proton recoil. The correction used in this work 
was evaluated numerically using the expression taken from a 
review article by Maximon [29]. 




where 



and E^ are the electron energies before and 
after scattering, 

a is the fine structure constant, 

m is the electron mass, 

2 . 

q IS the four momentum transfer squared, and 
L„ (x) = ~ f — — dt is the Spence Function. 

^ 'n t 



2 . Bremsstrahlung Correction 

When an electron passes through material, it will 
radiate energy through collisions with atomic electrons and 
nuclei other than the target nucleus of interest. This 
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energy loss which is referred to as bremss trahlung has been 
studied by Bethe and Heitler [30]. The energy loss per 
distance of material traversed is given by 



dE 

dx 



E X , 
o 



where is the radiation length, which is the mean distance 
over which the electron has its energy reduced by a factor 
of 1/e. 

Following Butcher and Messel [33] the radiation 
length for a multi-element target was calculated using 




where X^ is given by 



>=0 = 



A. 

1 



4a N 



r^ Z(Z+0.8) 



iln{183 Z 



where is Avogadro's number, 

r is the classical electron radius, 
e 

is the density in g/cm^ of the i^^ element, and 
A^ is in grams /mole of the i^^ element. 

Hofstadter obtained the bremsstrahlung correction 
factor by integrating the probability function P (a ) given 
by Bethe and Heitler. 



P(a) 



^(t/iln?) 
r (l+t/Jln2) 



where P(a) is the probability that an electron traversing 
t radiation lengths will have an energy loss less than a 
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times its initial energy (aside from the ionization loss and 
the scattering from the nucleus of interest) . The brems- 
strahlung correction factor is 



K 



b 



exp 



T 

X ln2 
o 



iln 




/ 



where T is the thickness of the target (this includes 

the two steel 1 mil windows as well as the gas) 
in g/cm2, and 

E is the geometric mean electron energy in the 
target. 

In this experiment the bremsstrahlung correction is 
about 2%. 

3. Ionization Correction 

Electrons also lose energy by excitation and ioniza- 
tion in the gas and in the gas target windows. Landau [34] 
calculated the average energy that a mono-energetic beam 
would lose, and the spread of the energy distribution of the 
beam after passing through a thin target. From-this distri- 
bution it is possible to calculate the fraction of electrons 
which are not detected because of energy degradation through 
ionization . 

The fraction, ij^(AE), of the beam one does not detect, 
because of ionization losses is 



E_ , E. 

r 2 / f 

tMAE) = J N(E) dE/J n(E) dE , 

o / o 



where E 2 is the energy corresponding to the energy of the 
last spectrometer setting. 



E^^^ is the incident electron energy, and 
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N(E) is the number of electrons detected at energy E. 

Landau [34] derived an approximation for ^(AE) , as 
have others [35], which is used in this work; 

ij;{AE) ^ 1/w , 



where the value of w as provided by Borsch-Supan [36] is 



0) Si 




S,nA +h\ 

X+1 j ' 



where A = Y~1 = -0.4228..., 

Y is Euler's constant, 

X = AE/C - 0.05 

I n . z . 

C = 0.154 ^ T MeV, and 

i i i 

2 

T = target thickness in g/cm , 

n • is the number of atoms of the i^^ element per 
molecule of the target, and and are 
the atomic number and atomic weight of the 
ith element. 

The ionization correction is thus given by 
= [1 - Ip (AE) ]"^ . 



For both the bremsstrahlung and ionization correc- 
tions, the target was treated as a homogeneous molecular 
compound of steel and either H 2 or D 2 , whereas in actuality 
the target consisted of H 2 or D 2 gas (2 inches thick) at 
77°K and 150 psig and two 1 mil steel windows. The molecu- 
lar formulas used to approximate the target were either 
Fe H 2 g or Fe D 2 g depending on the gas used. 
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D. ERROR ANALYSIS 



1 . Systematic Errors 

In this experiment the proton and deuteron cross 
sections were taken with the same geometry, current, and 
final electron energies. At most a time delay of twelve 
hours between the end of the hydrogen and the start of the 
deuterium run enabled the proton and deuteron data to be 
taken within a 48 hour time span. Consequently, in the 
ratio of the deuteron to proton cross sections, many of the 
systematic errors cancelled; e.g. solid angle, scattering 
angle, SEM efficiency, target thickness, and absolute 
spectrometer efficiency. 

A systematic error which would not cancel in this ratio 
experiment is that which arises from impurities in the 
hydrogen and deuterium gases. The gases used in this experi- 
ment were respectively 99.5% deuterium and 99.995% hydrogen, 
as determined by their supplier [37]. However, the equip- 
ment used to establish these purities was able to detect 
only 0.3% or more contamination of in H 2 or H 2 in D 2 • 

The data taken in this experiment used the same bottle of 
deuterium and two different bottles of hydrogen. Consequently, 
the hydrogen data may have varying amounts of possible con- 
taminants depending on which bottle was used. However, if 
there were a difference in the contamination of H 2 or D^i 

this difference would have caused a non zero intercept of 

2 

the plot of Gg versus q . This difference was not evident 
^ 2 

as the value of (q =0) is zero to within statistical 

n 

accuracy for all deuteron models tested. 
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2. Random Errors 



The statistical error is due to the Poisson distri- 
bution of the number of electrons scattered for a given 
number of incident electrons. This error was calculated by 
a subroutine in the BERHFORM and TRADFORM programs for each 
individual peak and included the statistical errors made in 
determining the background subtraction. 

If X is the cross section, the following may be 

written 



M 

X = (Constant) • Z (N, - B) , 

k=l 

where is the total counts (count rate and voltage 

corrected) in the k^h spectrometer setting, 

B is the average background to be subtracted and 
is assumed to be the same for each spec- 
trometer setting, and 

M is the number of spectrometer settings. 



It follows that if N = N. where N is the total number of 

k 

counts, then: 



X - (Constant) * (N - M-B) 



The variance of the cross section is 



2 2 
2 ^ ^ 2 ^ 2 

X 9N N 9B '^B ' 



0^ = (Constant)^ (N + M^0g) , 

2 

where is the variance in the average background, and is 
equal to B^. Let B^ = B^ (where B^^ is equal to the back- 
ground at the i'^*^ spectrometer setting) and let J be the 
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number of spectrometer settings above the peak at which 

J 

Z B . 

i_l 1 

background was taken. Since B = , then 




which in turn implies that 

0^ = (Constant) ^ (N + , 

A J 

and the fractional standard deviation of the cross section is 

2 

N + M B/J 
(N-MB)^ 

Besides the statistical (Poisson) counting errors, the 
random errors in this experiment originate from uncertainties 
in determining the density of gas in the target, variations 
in SEM efficiency during a run, energy shifts of the inci- 
dent beam, leakage currents in the beam current integrators, 
and small errors in the radiation corrections resulting from 
errors in determining AE . 

Uncertainties in determining the density of the gas in 
the target were caused by the variation of pressure in the 
target, that is, the pressure dropped while the liquid 
nitrogen was being pumped and then rose again and stayed 
steady until the next nitrogen filling. These variations 
were almost identical in the hydrogen and deuterium runs 
and were about 0.3% of the measured pressure. The average 
pressure for each peak in the hydrogen and deuterium runs 
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was determined separately by either taking the pressure at 
the start of the nitrogen fill or at the end of the fill 
depending on which pressure was more constant. If the 
pressure at the start (completion) of the nitrogen fill was 
used for hydrogen, then the pressure at the start (comple- 
tion) of the nitrogen fill was used for deuterium. The 
errors in the measured cross sections due to the uncertainty 
in determining the density of the gas in the target varied 
from run to run, but were always less than 0.1%. 

The efficiency of the SEM was measured four times during 
one data run on each gas, one on each of the four elastic 
peaks. On each peak, five or more measurements were made. 
The fractional standard deviation of a single measurement 
was usually less than 0.2%, giving a standard deviation of 
the mean of the SEM readings of less than 0.09% for the four 
readings . 

Due to variations in line voltage, the beam energy was 
not constant, and as a result the scattered energy also 
changed slightly from time to time throughout a run. Also 
there was a slight instability in the spectrometer power 
supply. The instabilities had the effect of shifting the 
counts from one channel to another on the counting ladder. 
This variation in the data was most noticeable on the steep 
slopes of the elastic peaks where a low count in one channel 
was accompanied by a high count in the adjacent channel. 

It was found that by closing the energy defining slits more, 
the variation in the counts due to shifts in energy was 
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almost negligible. At the same time, it was noticed that 
the more the energy defining slits were closed the higher 
the background became. Thus, the ratio of real counts to 
background counts became less, which in turn increased the 
errors on the measured cross section. The position for the 
energy defining slits was adjusted to give the minimum 
combined contribution to the error in the cross section 
from the effects of energy shifts and increased background. 
Even after this was done, in some cases there still existed 
a need to adjust for obvious energy shifts. This was 
accomplished by smoothing the data points to the general 
elastic peak line shape. It is estimated that in the 
worst case these energy shifts could contribute an uncer- 
tainty in the average measured cross section of as much as 
0.2%, which the usual error in the measured cross section 
due to these energy shifts was about 0.1%. 

In the circuits used to measure and integrate the beam 
current there existed slight leakage currents. If these 
leakage currents were not constant throughout a run, then 
they would introduce errors. These currents were measured 
and found to be constant as well as negligible [5]. As a 
result the only error attributable to these leakage currents 
existed when the electron beam current varied. The electron 
current was usually held constant to within 20%. This 20% 
variation would cause an error in the ratio of the measured 
cross sections of about 0.02%. 
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Another source of uncertainty is the inability to 
determine the exact position of the peak of the elastic 
spectrum from the plotted data. This uncertainty is 
usually about 0.01 MeV, and causes an error in the Schwinger, 
bremsstrahlung, and ionization corrections. The error in 
the measured cross section due to this AE uncertainty for 
all three radiation corrections was of the order of 0.05%. 

An uncertainty in the final energy at which one terminates 
the integration of the elastic spectrum, contributes approx- 
imately a 0.03% error to the measurement of the cross 
section . 

A list of all the errors discussed can be found in 

Table IV-1. Excluding statistical error (since it varies 

so much from run to run) , the average error in the cross 

section ratio, , is about 0.24%. The total 

Avg. Others 

error of the weighted mean ratio R was then formed from 



2 2 2 

Sf . , = S ^ 
total stat Av. Others 



(IV-11) 



One significant source of error has been reduced since 
Topping [5] performed his ratio experiment. It was found 
that if quadrupoles 3 and 4 were used instead of 1, 2, 3 and 
4 that the background was effectively 50% lower. The most 
plausible explanation for this reduction of background is 
that quadrupoles 1 and 2 defocused the electron beam, caus- 
ing the beam to hit the beam pipe and thereby increased the 
background. The data in this experiment was collected with 
only quadrupoles 3 and 4 on, and as a result, the peak to 
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TABLE IV-1 



RANDOM AND SYSTEMATIC ERRORS IN CROSS SECTION DATA 



Errors 


Average Values of Errors 


Statistical 


0.3% 


Gas Density 


0.1% 


SEM Efficiency 


0.09% 


Energy Shifts 


0. 1% 


Leakage Currents 


0.014% 


E Measurement 
(Radiation Corrections) 


0.05% 


E Measurement 
(Integration) 


0.03% 


Average Error 

(excluding counting statistics) 


0. 24.% 



background ratio of counts was much higher than in the 
previous work [5] . 

The uncertainty of Gg arises from the errors made in 

n 

measuring the individual cross sections. Let R = 

where Op and are the experimental cross sections for the 

proton and deuteron respectively. Let S and S, be the 

p Q 

standard deviation of -the proton and deuteron cross sections 
respectively. G„ is determined from 

rj 

_ n _ 

Gg (q^) , (IV-1 2) 

P 



n 



i^E j % 

V d py E 



- 1.0 



and G„ /G„ can be determined from 
^d ^p 
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Solving for the standard deviation in R, denoted by 
due to the errors in both cross sections, one obtains; 



2 2 
^R = ^ 




(IV-14) 



The standard deviation of the ratio G„ .G„ , denoted by 

Ij n hi 

d p 

S_ /G , is obtained from Equation IV-14, 

Gd P 



S,. /G = 

""d P 



^mott 


"Mr 




11 


mott 


M 


is (qp> . 


1 J 



R 



2 

From Equation IV-12, one can determine the standard 

deviation in G„ , denoted by S_ , to be: 

n 



(IV-15) 
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(q d) 



E 



E 



n 



'E 



E 



R 
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(IV-16) 



In order to combine the results obtained from four different 

channels in one run, let R. denote the ratio of cross sec- 

1 

th 

tions in the i^ channel and let S denote the standard 

th ^ 

deviation of the ratio in the i^ channel. It follows that 

the weighted mean ratio R can be calculated from 

m 



7 

2 

i=4 



R = _ 

m 7 



1 



2 

i=4 



(IV-17) 
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R. 
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where the standard deviation of this weighted mean ratio is 
calculated from 



S 



R 



m 




(IV-18) 



and are the experimental values of the ratio and 

m 

of the standard deviation to that ratio. 
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V. EXPERIMENTAL RESULTS 



A. NEUTRON CHARGE FORM FACTORS 



The 21 runs by Topping [5] and the twelve runs by the 
authors are reported here. Each run consisted of four 
measurements of R, the ratio of the elastic electron- 
deuteron cross section to the elastic electron-proton cross 
section. Of these 132 measurements of R, 9 were not 

included due to reasons given in Ref. [5]. 

2 

Table V-1 gives the q transferred to the deuteron, 
the scattering angle, the experimental ratio R^^, and the 
standard deviation of R determined both for the total 



m 

error and the statistical (Poisson) error for a given run. 
Also included are quantities that were determined by cal- 
culations discussed in Section II. The ratio of the 
deuteron electric form factor to the proton electric form 
factor, derived using relativistic corrections by both 
Friar and Gross, is also listed in Table V-1. 

Tables V-2,3,4 and 5 present the neutron charge form 
factors obtained using the deuteron models LF1,5, and 15 
with the relativistic corrections of Gross and Friar. The 
deuteron structure factors used are presented in Appendix D. 
The best estimates of the error in each of these 33 values 

of G„ for each LF model and relativistic correction is the 
E 

n 

error obtained using from Equation IV-11. An average 

2 

G„ and the associated S. . ^ at each q were calculated 

E total ^ 

n 

using Equations V-1 and V-2 respectively. 
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i (Sf^.) 



tot i-i 



-U 



(V-2) 



A plot of the average G^. (with its versus q is pro- 



n 



vided for each of three LF models and both relativistic 
corrections . 

All values of G„ were calculated on a Wang Model 700-C 

n 

calculator using programs written by the authors. The best 



2 2 

estimate of the slope of G„ versus q , d/dq G„ , was 

n 



n 



obtained from a linear fit weighted by for all values 



of G^ 



For comparison purposes the linear fit was also 



n 



forced through zero for each plot. 

2 

The average value of G„ at each q , using Lomon- 

£j 

n 

Feshbach Model #1 is given in Table V-2. 

2 

The average value of G„ at each q , using Lomon- 

Feshbach Model #5 is given in Table V-3. 

2 

The average value of G at each q , using Lomon- 

lit 

n 

Feshbach Model #15 is given in Table V-4. 



The value of G„ calculated from Lomon-Feshbach Model 

hj 

n 

#15 and its Sq for each of the 33 runs is given in 
^n 

Table V-5. 
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B. ROSENBLUTH PLOT 
2 -2 

At q = 0 . 4 f ra deuteron and hydrogen data v/ere taken 

at 90° and 120° so that a Rosenbluth plot could be con- 
structed. The theoretical interpretation of the Rosenbluth 
plot is discussed in Appendix A. Table V-6 gives the 
experimental data and derived quantities used in the Rosen- 
bluth plot given in Figure V-7. The use of the scaling 
law is compatible with the data. 



C. CHARGE RADIUS OF THE DEUTERON 

Using Method II given by Schumacher and Bethe [7] , the 

deuteron structure rms radius was calculated. A more 

detailed explanation is given in Appendix B. The structure 

radius of the deuteron for LFl, LF5, and LF15 was provided 

by Schumacher and Bethe [7] and was used in a plot of the 

2 



experimental 



, 2 ^E 
dq n.. 



(deuteron structure 

q^=0 ""d 

radius) . This plot is Figure V-8 and as predicted by 



versus r„ 
,2=r E. 



Schumacher such a plot gives a straight line for all models 

with the correct binding energy. The intersection of this 

d 



line with the value of 



L , 2 ^E , 

Ldq n J 



q =0 



determined by the very 



accurate neutron-electron scattering work yields r^ . The 

^d 

value for the deuteron structure rms radius was determined 
from Figure V-8 to be r„ = 1.9665 ± 0.0045. 

X!j 
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TABLE V-1 



EXPERIMENTAL RATIOS AND RESULTS OF 
MODEL INDEPENDENT CALCULATIONS 



2 


6 


Exp. Ratio 
d 

R = ^ 

P 

exp 


Std. Dev. 
of R 


aP 

mott 




(qp) 


P 




Total and 
(stat . ) 


"^mott 




0.05 


60° 


0.9654 


0.0046 




0.30269 


0.99400 


0.99404 








(0.0040) 


0.31356 










75° 


0.9724 


0.0038 




0.16837 


0.99396 


0.99404 








(0.0030) 


0.17578 








0.10 


75° 


0.9390 


0.0042 




0.81525 


0.98789 


0.98815 








(0.0035) 


0.86757 










75° 


0.9331 


0.0034 




0.81525 


0.98789 


0.98815 








(0.0024) 


0.86757 










75° 


0.9262 


0.0046 




0.81525 


0.98789 


0.98815 








(0.0038) 


0.86757 








2 




(Gross) 




(Gross ) 

G_ Total 


1 


(Friar) 

Total 




0 


20 


rP 

M 










(Friar) 


G„ i(stat) 


i 


— (stat) 






r 






p Error 


3 Error 


0.05 


60° 


1.00001 


0.99338 


0. 


.9685±0.0023 


0.968510.0023 






1.00001 






(0.0020) 




(0.0020) 




75° 


0.99997 


0.99114 


0. 


.969310.0019 


0.969310.0019 






0.99996 






(0.0015) 




(0.0015) 


0.10 


75° 


0.99993 


0.98230 


0. 


.947510.0021 


0.947510.0021 






0.99992 






(0.0018) 




(0.0018) 




75° 


0.99993 


0.98230 


0. 


.944510.0017 


0.944510.0017 






0.99992 






(0.0012) 




(0.0012) 




75° 


0.99993 


0.98230 


0, 


.941010.0023 


0.941010.0023 






0.99992 






(0.0019) 




(0.0019) 
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TABLE V-1 (Continued) 



2 


e 


Exp. Ratio 

d 

o 

R = ^ 

0P 

exp 


Std. Dev. 
of R 


mott 


P 


S «Jd> 

P 




Total and 
(Stat. ) 


mott 


0.10 


75° 


0.9221 


0.0040 




0.81525 


0.98789 


0.98815 








(0.0031) 


0.86757 








0.20 


60° 


0.8204 


0.0051 




0.70815 


0 . 97595 


0.97653 








(0 . 0045) 


0.76120 










60° 




0.0076 




0.70815 


0.97595 


0. Q7653 








(0.0072 


0.76120 










75° 


0.8402 


0.0049 




0.30041 


0.97582 


0.97653 








(0.0042) 


0.42650 










90° 


0.8323 


0.0049 




0.22434 


0.97569 


0 . 97653 








(0.0042) 


0.24871 










90° 


0 . 8365 


0.0052 




0.22434 


0 . 97569 


0.97653 








(0.0046) 


0.24871 








2 

‘Jd 


0 


(Gross) 


rP 

M 


s 


(Gross ) 

, Total 

'd ^ and 

■ (Stat) 


(Friar) 

G Total 

^ rs ^ 






(Friar) 

Mr 






■ (Stat) 










'p Error 


P 


Error 


0.10 


75° 


0.99993 


0.98230 


0. 


9389±0 . 0020 


0.938910.0020 






0.99992 






(0.0016) 




(0.0016) 


0.20 


60° 


1.00001 


0.97359 


0. 


884810.0028 


0.884810.0028 






1.00001 






(0.0024) 




(0.0024) 




60° 


1.00001 


0.97359 


0 . 


888710 . 0041 


0 . 888710 . 0041 






1.00001 






(0.0039) 




(0.0039) 




75° 


0.99985 


0.96493 


0 . 


892110.0026 


0.892110 .0026 






0.99985 






(0.0022) 




(0.0022) 




90° 


0.99896 


0.95129 


0 . 


887110.0026 


0.887110.0026 






0.99896 






(0.0022) 




(0.0022) 




90° 


0.99896 


0 . 95129 


0 . 


889310.0028 


0 . 889310 . 0028 






0.99896 






(0 . 0024) 




(0.0024) 
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TABLE V-1 (Continued) 



2 


6 


Exp. Ratio 
d 

0 

R = 

aP 

exp 


Std. Dev. 
of R 


cP 

mott 


p 1 




=E <‘3d> 

p 




Total and 
(Stat. ) 


d 

^mott 


^E 

P 


0.20 


120° 


0.8068 


0.0118 




0.72058 


0.97555 


0.97653 








(0.0106) 


0 . 81829 










120° 


0.8280 


0 . 0102 




0.71893 


0.97550 


0.97653 








(0.0099) 


0.81637' 










120° 


0.8057 


0.0082 




0 . 71893 


0.97550 


0.97653 








(0.0079) 


0.81637 








0.25 


75° 


0.7987 


0.0056 




0.30690 


0.96982 


0 .97082 








(0.0051) 


0 . 33890 








0.30 


75° 


0.6973 


0.0049 




0.27183 


0.96625 


0.96516 








(0.0043) 


0.28073 










90° 


0.7718 


0.0057 




0 . 14386 


0 . 96365 


0.96516 








(0.0052) 


0.16338 








2 

‘Jd 


6 


(Gross) 

r 


rP 

M 




(Gross ) 

Total 

'd ^ and 

■ (Stat) 


(Friar) 

G Total 






(Friar) 






~ (Stat) 






Mr 






‘p Error 


p 


Error 


0.20 


120° 


0.99466 


0.89020 


0 


. 8901±0.0065 


0. 8901±0.0065 






0.99466 






(0.0058) 




(0.0058) 




120° 


0.99465 


0.88999 


0 


.9018±0.0055 


0 . 901810.0055 






0 . 99464 






(0.0054) 




(0.0054) 




12 0° 


0.99465 


0.88999 


0 


. 8896±0. 0045 


0. 889510.0045 






0.99464 






(0.0044) 




(0.0044) 


0.25 


75° 


0.99982 


0.95641 


0 


.8687±0.0030 


0.868710.0030 






0.99981 






(0.0028) 




(0.0028) 


0.30 


75° 


0.99978 


0.95135 


0 


.8426±0.0030 


0.842610.0030 






0.99977 






(0 . 0026) 




(0.0026) 




90° 


0.99845 


0.92817 


0 


.8538±0.0031 


0.853810.0031 






0.99844 






(0.0029) 




(0.0029) 
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TABLE V-1 (Continued) 



2 


e 


Exp. Ratio 

d 

a 

R = -^ 

0P 

exp 


Std. Dev. 
of R 


*^mott 


P 




S '-Sd' 

p 




Total and 
(Stat. ) 


_d 

mott 


0.30 


90° 


0.7529 


0.0060 




0.14386 


0.96365 


0.96516 








(0.0055) 


0.16338 










90° 


0.7558 


0.0070 




0.14386 


0.96365 


0.96516 








(0.0066) 


0 . 16338' 








0.35 


90° 


0.7192 


0.0054 




0 . 12122 


0.95767 


0.95955 








(0.0048) 


0.13916 










90° 


0.6596 


0.0056 




0.13341 


0.96112 


0 . 95955 








(0.0051 




0. 13<316 








0.40 


90° 


0.6815 


0.0056 




1.04376 


0.951718 


0.95400 








(0.0050 




1. 210236 










90° 


0.6829 


0.0058 




1.04376 


0.95172 


0.95400 








(0.0053) 


1.21024 










e 


(Gross) 

r 


< 




(Gross ) 

. Total 

^ 


(Friar) 

G Total 


(fm 




(Friar) 

Mr 


G„ " (Stat) 




■ (Stat) 










‘p Error 


P 


Error 


0.30 


90° 


0.99845 


0.92817 


0 


. 8433±0.0034 


0.843310.0034 






0.99844 






(0.0031) 




(0.0031) 




90° 


0 .99845 


0 . <32817 


0 


. 844°±0.0039 


0. 844910.0039 






0.99844 






(0.0037) 




(0.0037) 


0.35 


90° 


0.99819 


0.91695 


0 


. 824210.0031 


0.824210.0031 






0.99818 






(0.0027) 




(0.0027) 




90° 


0.99819 


0.92340 


0 


. 8281 + 0.0035 


0.828110.0035 






0.99818 






(0.0032) 




(0.0032) 


0.40 


90° 


0.997933 


0 . 90 5924 


0 


. 802710 . 1 


0033 


0.802710.0033 






0.997918 






(0.0030) 




(0.0030) 




90° 


0.997933 


0.905^24 


0 


.803610. i 


0034 


0.803610.0034 






0 .997918 






(0.0031) 




(0.0031) 
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TABLE V-1 (Continued) 



2 


6 


Exp. Ratio 

d 

a 

R = 

qP 

exp 


Std. Dev. 
of R 


*^mott 






=E <-5a> 

P 




Total and 
(Stat. ) 


mott 


^E ' 
P 


0.40 


120° 


0.6430 


0 . 0062 




0.32879 


0.95117 


0.95400 








(0.0057 


) 


0.39489 










12 0° 


0.6266 


0.0061 




0.32879 


0.95117 


0.95400 








(0.0056 


) 


0.39489 










120° 


0 .6300 


0.0062 




0.32879 


0.95117 


0.95400 








(0.0057 


) 


0.39489 








0.477 


90° 


0.6324 


0.0061 




0.85522 


0.94262 


0.94550 








(0.0056 


) 


1.00679 










90° 


0.6389 


0.0051 




0.85522 


0.94262 


0.94550 








(0.0045 


) 


1.00679 










90° 


0.6330 


0.0060 




0 . 85522 


0.94022 


0.94500 








(0.0055 


) 


1.00679 








2 

‘^d 


e 


(Gross) 


rP 

M 


S 


(Gross ) 

. Total 

[d and 

, “ (Stat) 


(Friar) 

G Total 






(Friar) 

Mr 




P 


“ (Stat) 










'p Error 


Error 


0.40 


120° 


0.98938 


0.79894 


0, 


. 8118±0.0039 


0.811810.0039 






0.98986 






(0.0036) 




(0.0036) 




120° 


0.98938 


0.79894 


0 . 


.801410.0039 


0.801410.0039 






0.98936 






(0.0036) 




(0.0036) 




120° 


0 .98938 


0.79894 


0. 


.803610.0039 


0.803610.0039 






0.98936 






(0.0036) 




(0 . 0036) 


0.477 


, 90° 


0.99754 


0.88935 


0, 


.773910.0037 


0.773910.0037 






0.99752 












(0.0034) 




! 90° 


0.99754 


0.88935 


0, 


.777810.0031 


0.777810.0031 






0.99752 






(0.0027) 




(0.0027) 




90° 


0.99754 


0.88935 


0 


.774210.0037 


0.774210.0037 












(0.0034) 




(0 .0034) 
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TABLE V-1 (Continued) 



2 


0 


Exp. Ratio 
d 

0 

R = -^ 

aP 

exp 


Std. Dev. 
of R 


*^mott 


p ^ 






Total and 
(Stat. ) 


*^mott 


0.50 


105° 


0 . 6033 


0.0057 

(0.0052) 


0.46122 

0.55662 


0.93952 


0.94307 




105° 


0 . 6042 


0 . 0057 
(0.0052) 


0.46122 

0.55662 


0.93952 


0.94307 




105° 


0 .6044 


0 .0058 
(0.0053) 


0.46122 

0.55662 


0.93952 


0.94307 




10 5° 


0.6015 


0.0060 

(0.0055) 


0.46122 

0.55662 


0.93952 


0.94307 


2 

‘^d 


0 


(Gross) 


rP 

M 


s 


(Gross ) 

Total 

'd and 

■ (Stat) 
'p Error 


(Friar) 

G Total 

^ ^ 






(Friar) 

Mr 


s 


^E 

P 


■ (Stat) 
Error 


0.50 


105° 


0.99368 

0.99366 


0.83642 


0 


. 7677±0. 0036 
(0.0033) 


0. 7677±0.0036 
(0.0033) 




105° 


0.99368 

0.99366 


0.83642 


0 


. 7683±0 . 0036 
(0.0033) 


0 . 7683±0. 0036 
(0.0033) 




105° 


0.99368 
0 . 99366 


0.83642 


0 


.7685±0.0037 
(0.003f ) 


0.7685±0.0037 

(0.0034) 




105° 


0.99368 

0.99366 


0 . 83642 


0 


. 7666±0. 0038 
(0.0035) 


0.7666±0.0038 

(0.0035) 
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TABLE V-2 



AVERAGE G_ AT 
E 

n 


EACH q^ USING 


LOMON-FESHBACH MODEL #1 


2 

q 


Gross' Relativistic 
Correction 




Friar's Relativistic 
Correction 




n 


En 








Sr 


0.05 


0.0002 


0.0015 




0.0002 




0.0015 


0.10 


0.0037 


0.0011 




0.0037 




0.0011 


0.20 


0.0027 


0.0013 




0.0027 




0.0013 


0.25 


0.0070 


0.0034 




0.0069 




0.0034 


0.30 


0.0076 


0.0019 




0.0075 




0.0019 


0.35 


0.0093 


0.0027 




0.0092 




0.0027 


0.40 


0.0084 


0.0020 




0.0081 




0.0020 


0.477 


0.0098 


0.0025 




0.0095 




0.0025 


0.50 


0.0106 


0.0023 




0.0102 




0.0023 


A linear 


fit to the 


data in Table 


V-2 


gives : 






S = 


(0.0003 ± 


0.0010) + (0. 


0208 


+ 0.0040)q^ 




n 


for Gross' 


correction and 








^E = 


(0.0004 + 


0.0010) + (0. 


0200 


± 0.0040)q^ 




n 


for Friar' 


s correction. 










A linear 


fit of the 


same data forced 


through 


zero 


gives : 


Gj, (q 
n 


= (0.0218 ± 0.0022)q^ 


for 


Gross ' 


correction and 


n 


= (0.0213 ± 0.0022)q^ 


for 


Friar ' s 


correction. 
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0.02-DERlVED FROM LF1 AND 
GROSS' RELATIVISTIC CORR. 
Q_ vs. 




o 

a 



o 

d 



o 



LjJ 
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Figure V-1 
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TABLE V-3 



AVERAGE 


AT EACH 

Ij 

n 


USING 


LOMON-FESHBACH MODEL #5 


2 

q 


Gross' Relativistic 
Correction 




Friar's Relativistic 
Correction 




n 


^Ge 

^n 




n 


^n 


0.05 


0.0002 


0.0015 




0.0002 


0.0015 


0.10 


0.0037 


0.0011 




0.0037 


0.0011 


0.20 


0.0028 


0.0013 




0.0027 


0.0013 


0.25 


0.0070 


0.0034 




0.0069 


0.0034 


0.30 


0.0076 


0.0019 




0.0075 


0.0019 


0.35 


0.0094 


0.0027 




0.0092 


0.0027 


0.40 


0.0084 


0.0020 




0.0082 


0.0020 


0.477 


0.0099 


0.0025 




0.0096 


0.0025 


0.50 


0.0107 


0.0023 




0.0103 


0.0023 


A linear fit 


to the data 


in Table 


V-3 


gives : 




S (^3") = 

n 


= (0.0003 ± 
for Gross' 


0.0010) + (0.0210 ± 0 
correction and 


.0040)q^ 


Ge (q2) = 
n 


= (0.0004 ± 
for Friar' 


0.0010) + (0. 
s correction. 


0202 ± 0 


. 0040)q^ 


A linear fit 


of the same 


data forced 


through 


zero gives: 


S ^^3') = 

n 


: (0.0220 ± 


0.0022)q^ 


for 


Gross ' 


correction and 


Ge (q") = 
n 


: (0.0214 ± 


0.0022)q^ 


for 


Friar ' s 


correction. 
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Figure V-2 



TABLE V-4 



AVERAGE Gg AT. EACH USING LOMON-FESHBACH MODEL #15 
n 



2 

q 


Gross' Relativistic 
Correction 




Friar's Relativistic 
Correction 


n 


^n 


^E 

n 


Sp 

^En 


0.05 


0.0000 


0.0015 




-0.0000 


0.0015 


0.10 


0.0033 


0.0011 




0.0033 


0.0011 


0.20 


0.0020 


0.0013 




0.0020 


0.0013 


0.25 ' 


0.0061 


0.0034 




0.0060 


0.0034 


0.30 


0.0066 


0.0019 




0.0066 


0.0019 


0.35 


0.0082 


0.0027 




0.0080 


0.0027 


0.40 


0.0070 


0.0020 




0.0068 


0.0020 


0.477 


0.0083 


0.0025 




0.0080 


0.0025 


0.50 


0.0090 


0.0023 




0.0086 


0.0023 


A linear fit 


to the 


data in Table 


V-4 


gives : 




,G^ (q2) . 


= (0.0002 ± 0.0010) + 


(0. 


0177 ± 0.0040)q^ 


n 


for Gross ' correction and 




Ge (q2) ^ 


= (0.0003 ± 0.0010) + 


(0. 


0170 ± 0.0040)q^ 


n 


for Friar's correction. 






A linear fit 


of the 


same data forced 


through zero 


gives : 


G (q2) 


= (0.0185 ± 0.0022)q^ 


for 


Gross' correction and 


i-i 

n 












Ge (q^) 


= (0.0180) ± 0.0022)q^ 


) 

for Friar's correction. 


n 
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Figure V-3 




65 



Figure V-4 
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Figure V-6 




f 




TABLE V-5 



AND ^ FOR EACH OF THE 30 RUNS 
E tot 

n 

CALCULATED FROM LOMON-FESHBACH MODEL #15 



2 

q 

(fm”^) 


Gross' Relativistic 
Correction 


Friar's Relativistic 
Correction 


^E 

n 




^E 

n 


^n 


0.05 


0.0003 


0.0019 


0.0003 


0.0019 




-0. 0005 


0.0024 


-0.0005 


0.0024 


0.10 


0.0078 


0.0022 


0.0078 


0.0022 




0.0010 


0.0025 


0.0010 


0.0025 




0.0047 


0.0018 


0.0047 


0.0018 




-0.0012 


0.0021 


-0.0012 


0.0021 


0.20 


-0.0027 


0.0030 


-0.0028 


0.0030 




0.0015 


0.0045 


0.0014 


0.0045 




0.0052 


0.0029 


0 .0052 


0.0029 




-0.0002 


0.0029 


-0.0003 


0.0029 




0.0021 


0.0030 


0.0021 


0.0030 




0.0030 


0.0072 


0.0030 


0.0072 




0.0159 


0.0061 


0.0158 


0.0061 




0.0025 


0.0050 


0.0024 


0.0050 


0.25 


0.0061 


0.0034 


0.0060 


0.0034 


0.30 


0.0153 


0.0036 


0.0151 


0.0036 




0.0032 


0.0039 


0.0030 


0.0039 




0.0050 


0.0045 


0.0049 


0.0045 




0.0024 


0.0034 


0.0023 


0.0034 


0.35 


0.0061 


0.0036 


' 0.0060 


0.0036 




0.0108 


0.0041 


0.0106 


0.0041 


0.40 


0.0050 


0.0039 


0.0048 


0.0039 




0.0060 


0.0041 


0.0058 


0.0041 




0.0159 


0.0047 


0.0157 


0.0047 




0.0034 


0.0046 


0.0032 


0.0046 




0.0060 


0.0047 


0.0058 


0.0047 


0.477 


0.0061 


0.0046 


0.0058 


0.0046 




0.0110 


0.0038 


0.0107 


0.0038 




0.0066 


0.0045 


0.0063 


0.0045 


0.50 


0.0089 


0.0045 


0.0086 


0.0045 




0.0096 


0.0045 


0.0092 


0.0045 




0.0098 


0.0046 


0.0095 


0.0046 




0.0074 


0.0047 


0.0071 


0.0047 
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The slope from the scaling law is 0.00097. 
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Figure V-7 




8 



Figure V-8 
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VI. CONCLUSIONS 



As a result of the 123 measurements of the ratio of the 

elastic electron-deuteron cross section to the elastic 

electron-proton cross section, nine values of at dif- 

2 

ferent values of q were determined for each deuteron wave 

2 

function considered. From these values of G„ and q , a 

hi 

n 



value of 



\—2 S 1 

Ldq^ ni 



q ^=0 



was determined. 



Relativistic corrections proposed by both Friar [8] and 

Gross [6] were investigated. It was found that, for a given 

d 



deuteron wave function. 



, 2 

dq n ' 



was smaller using the 



q =0 

Friar correction. Although this pattern was consistent, it 
was also found that the difference in the two slopes was 
always well within one standard deviation. No attempt was 
made to select one form of the relativistic correction as 
best because the difference was so small. 

All wave functions considered (LFl, 5, 15) were in 
agreement with the very accurate thermal neutron slope of 



(0.0193 ± .0004) fm . The 



values of 2 .. obtained 

Ldq2 "Jq =0 



using Gross' relativistic corrections were: for LFl (0.0218 ± 

± 0.0022), for LF5 (0.0220 ± 0.0022) and for LF15 (0.0185 
± 0.0022). The values obtained using Friar's corrections 
were: for LFl (0.0213 ± 0.0022), for LF5 (0.0214 ± 0.0022), 

and for LF15 (0.0180 ± 0.0022). LF15 was found to be in 
closest agreement with the thermal neutron slope. 
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A value for the structure radius of the deuteron was 
also obtained from this data. This value was determined 
by comparing the slopes given by each wave function. The 
structure radius was determined to be r_ = (1.9665 ± 0.0045)fm. 

tj 

2 -2 

A Rosenbluth plot was made at q = 0 . 4 f m . The results 

of this plot are discussed in Appendix A and show agreement 

with the use of the conventional deuteron scaling law. 

If this experiment were to be repeated or extended it is 

2 

recommended that data at higher q be taken in order to 

4 

determine more accurately the coefficient of the q term in 
the expansion of G„ . In Appendix C a coefficient for the 

rj 

4 ^ 

q term was taken from work done at DESY [38] and it was 

2 -2 

found that at q =0.5 fm this term decreases G„ by 

Xj 

n 

about 10%. The effect on the slope of the G^ curve is also 

n 

shown to be about a 7% increase. 
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APPENDIX A 



THE ROSENBLUTH PLOT 



Based on the one photon exchange, the cross section for 
electron-proton scattering was derived by Rosenbluth [10]. 
Using the same assumption (one photon exchange) Jankus [14] 
calculated the deuteron cross section. Jankus' result is 
given in Equation II-7 and is reproduced below 



da 

dJ^ 



da 



dfi 



fmott 



_2 ^ 8 ^2 _2 , 2 _2 

^E,"^ 9 ^0, 3 Si, 

d ' d d 



1 + 2 ( 1+n ) tan^ -|- 



where G, can be written as 
d 



- 

^d “ 



’2 , 2 . , 8 ^2 _2 ^ 2 ^ „2 

"E, ^ 9 ^ ^Q, ■ 3 ^ 

d '^d d 






2 bI 

1+ 2 (1+n) tan |-J 



Since the deuteron cross section can be measured at different 

scattering angles while at the same four momentum transfer, 

2 2 0 

one can see that a plot of versus (1 + 2 (1+n) tan 2 ) 

. ■ 2 2 
should give a straight line with -j n as the slope and 

2 8 2 2 ^ 

+ E’ 0 G^ as the intercept. By plotting the deuteron 

^d ^ ^d 

^ ^ 2 , 2 6 
total form factor G^ versus [l+2(l+n)tan ^] , the adequacy 

of the conventional scaling law can be tested. Such a plot 

2 -2 

has been made for q = 0 . 4 fm and it can be found in 
Section V. 

One cannot use Equation II-7 directly since it contains 
the absolute cross section of the deuteron. However the 
deuteron cross section can be obtained by using the hydrogen 
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i 

i 



I 



\ 

( 

i 



i 

1 



i 



I 

( 



ii 



cross section along with the deVries b' fit. The deVries b' 
cit [19] is used because it is the best fit to the absolute 
proton cross section. It is used as a normalization factor 
with the experimental cross section obtained from the data. 
That is. 



da^ 



= N 



(deVries 



da^ 



fexp 



Therefore , 




N is a normalization factor for a given experiment at a 
given angle and energy. Since both the proton cross section 
and the deuteron cross section were measured under the same 
experimental conditions, the deuteron cross section is 



da 



d^ 



= N 



da 






(exp 



This implies that 

, , d 
da 



,dQ 



'doP' 






deVries 



da^ 



dfi 



/exp 



da 

dil 






(exp 



Now putting this in terms of the experimental ratio R, 
where R is given by 



R = 



'da' 

,da 



fexpi 



da^ 

d^l 



fexp 



one obtains 
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do 

1 



= R 



do^ 

dfi 



deVries 



and recalling that 

j \ 



do' 

dfi 



do 



dQ I d ' 

fmott 



one gets the final expression for G, as 






do^ 

dfi 



'deVries 



do 



dfi 



'mott 



R , 

m 



where 



^d 
doP 

j was calculated from the BERROTON program 

'deVries 



and 



do 






was calculated from the BERDEUTR program and 



mott 



the ratio R^^ is the weighted mean ratio given by Equation 
IV-17. 

2 

The statistical error in was calculated using the 
following equation; 



do* 



, dQ' /deVries 



do' 



dJ2 



R ' 

m 



mott 



where S is given by Equation IV-18. 

K. 

m 
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APPENDIX B 



RMS STRUCTURE RADIUS OF DEUTERON 

By substituting the deuteron charge distribution for the 
neutron charge distribution, it can be shown that an expres- 
sion similar to Equation 11-15 for the deuteron is 
applicable . 



2 4 

Gd(q ) = ~ 5 ^ 120 + ••• 



2 . 



(B-1) 



where r^ is the mean square radius of the deuteron charge 

distribution. In principle, the deuteron form factor could 

2 

be fit to this polynomial form and the slope of G , (q ) ^ 

<r 

2 2 d 

versus q in the limit as q approaches zero v;ould be - — — 



However, for the deuteron, the coefficient of the q term is 

4 

not well determined and the choice of greatly affects 

2 4 

the value of <r^>. In order to determine <r->,. data at 

d d 

higher momentum transfer is necessary and its determination 
requires knowledge of the coefficient of the q^ term. 

Fortunately, Schumacher and Bethe [7] have suggested a 

procedure to avoid this problem. The method consists of 

. 2 
using the experimental data and the function C„ (q ) associ- 

hi 

ated with a deuteron wave function to calculate A in the 
expansion 



(q ) = K 



B 



(B-2) 



n 



77 



If the deuteron structure factor is also expanded for 

£j 

2 

small q , the expansion is: 

Cj,(g2) = 1 - I <r^> q2 + <r^> q" + (B-3) 

2 

where <r„> may be defined as the mean square radius of the 

hi 

deuteron structure factor. 

Schumacher and Bethe observed that if the same data set 

is used to calculate A from various C„ functions, which have 

2 

different values of <r„>, the plot of A as a function of 

£j 

2 

<r> forms a straight line, if one uses only C„ functions 

li tj 

corresponding to deuteron models which have the correct 

binding energy. Then, if one requires that A have the value 

deduced from the neutron-electron scattering experiment, the 
2 

value of niay be determined, as shown in Figure V-8. 

2 

The error in <r > is fixed by the intersection of the pro- 

Ij 

jection of the errors of A, and the slope obtained using 

2 

models considered in this experiment. The value of 
is given in Section V. 

Finally, one may note that the mean square charge radius, 

2 

<r^> may be found from; 
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APPENDIX C 



DEPENDENCE OF 

n 

In Chapter IV, the slope of G„ was determined assuming 

hi 

2 

a linear fit of G versus q , but there was some evidence 

ill 

n 

that the inclusion of a quadratic term may be more appro- 
priate. Therefore, a fit with both linear and quadratic 
terms was made to the data, but it was discovered that 
such a fit produced a quadratic coefficient of little 
confidence . 

Since the value of the coefficient for the linear term 
is not strongly dependent on the value of the quadratic 
coefficient, for this data, a value for the quadratic 
coefficient was assumed using the work of Galster et al . 

4 

[39]. Galster al . obtained the value B = -0.0036 fm 
for the quadratic coefficient and this value was used in 
application to this work's data. A recalculation for the 
quadratic coefficient of the analytic expression, 

2 4 

Gg = Aq^ + Bq^ , 
n 

shows that the quadratic term contributes 5% to G„ at 

tj 

2-2 2-2 ^ 
q = 0.35 fm and 10% at q = 0 . 5 f m . With the assumed 

quadratic coefficient, the slope of G„ increased by 7% 

tj 

n 

from that obtained without the inclusion of the quadratic 
term. Also, the deuteron structure radius is decreased by 
0.15%. Note that the changes in the slope and the deuteron 
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structure radius are within the statistical accuracy of 
the value previously reported in Chapter IV. A comparative 
table of values utilized in the linear and quadratic fits 
can be found in Tables C-1 and C-2, 
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TABLE C-1 



COMPARISON OF AND 



2 

q 


LF5 




LF15 


(fm 


S S 

n n 


n 


S -Bq' 

n 


0.05 


0.0002 0.0002 


0.0000 


0.0000 


0.10 


0.0037 0.0037 


0.0033 


0.0033 


0.20 


0.0028 0.0029 


0.0020 


0.0021 


0.25 


0.0070 0.0072 


0.0061 


0.0063 


0.30 


0.0076 0.0079 


0.0066 


0.0069 


0.35 


0.0094 0.0098 


0.0082 


0.0086 


0.40 


0.0084 0.0089 


0.0070 


0.0076 


0.477 


0.0099 0.0107 


0.0083 


0.0091 


0.50 


0.0107 0.0116 


0.0090 


0.0098 


For LFl 


4 

assuming B = - 0.0036 fm 

A = 0.0232 ± 0.0022 






For LF5 


4 

assuming B = - 0.0036 fm 

A = 0.0232 ± 0.0022 






For LF15 


assuming B = - 0.0036 fm^ 

A = 0.0198 ± 0.0023 


- 






TABLE C-2 






COMPARISON OF THE DEUTERON STRUCTURE RADIUS OBTAINED 

FROM THE LINEAR AND QUADRATIC FITS TO 

G^ , ASSUMING GROSS RELATIVISTIC CORRECTION 
E 

n 




G = Aq2 

n 


Gp = Aq 
^n 


2 ^ u 4 
: + bq 




3.867 ± 0.018 


3.858 ± 


0.018 


r„ 1.9665 ± 0.0045 

^rms 


1.9642 ± 


0.0045 
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APPENDIX D 



COMPARISON OF GROSS' AND FRIAR'S RELATIVISTIC CORRECTION 



Both the Gross and Friar relativistic corrections 
multiply the right hand side of Equation II-7. The Gross 
relativistic correction factor is given by: 



1 




2 4 

_ _ 3 _ + 

2 4 ' 

4M 64 M 
P P 



(D-1) 



and Friar's correction factor is 




Expanding Friar's correction in a binomial series to second 
2 

order in ^ , the following is obtained: 

4M 

P 




The difference between the two relativistic correction 
factors is approximately given by the difference between 
Equation D-2 and D-1. If this difference is given the 
symbol A then. 



A 




64 



_3_ 

64 




(D-3) 



At 



2 

The difference is seen to increase with increasing q . 

2-2 2 
q = 0.5 fm , which is the highest q investigated, A 
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is 2.3 X 10 ^ with Friar's correction always being greater. 

2 -2 

Since the corrections are both about 0.997 for q = 0.5 fm , 
then the difference in these two corrections is about two 

5 

parts in 10 . Table D-1 gives a comparison of both correc- 

2 

tions for all q 's investigated 

Friar [8] also shows that C„, the electric structure 

Ij 

2 

factor of the deuteron, should not be evaluated at q but 

2 ' 

at an effective q given by 



q 



2 ' 




1 + 




(D-4) 



This also creates a very small change as seen in Tables D-2, 
D-3 and D-4. Friar's first relativistic correction opposes 
the tendency of the use of the effective momentum transfer, 
so the net difference between the use of the Friar and Gross 

4 

corrections is smaller than four parts in 10 . Also 

included in these tables is a comparison of the slope of 
2 

Gg versus q using each correction. Again the difference 
is seen to be negligible for this experiment. 
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TABLE D-1 



COMPARISON OF RELATIVISTIC CORRECTION FACTORS 



q" 

fm-2 


Friar ' s 
Relativistic 
Correction 


Gross ' 

Relativistic 

Correction 


^ Friar 
Gross 


0.05 


. 999972 


. 999972 


-.000000 


0.10 


.999447 


.999448 


-.000001 


0.20 


.998894 


.998896 


-.000002 


0.25 


.998618 


.998621 


-.000003 


0.30 


.998341 


. 998345 


-.000004 


0.35 


.998065 


. 998070 


-.000006 


0.40 


.997788 


.997796 


-.000007 


0.477 


.997360 


.997370 


-.000010 


0.50 


.997235 


. 997247 


-.000011 



TABLE D-2 

2 2 ' 

COMPARISON OF (q ) AND (q ) FOR LFl 


2 

q 


Cg(q') 


2 

q 


CE(q"') 




Cg(q') 


0.05 


.96882 


0.0500 


.96882 


-.00000 


0.10 


.93967 


0.1000 


.93969 


-.00001 


0.20 


. 88668 


0.1999 


.88673 


-.00006 


0.25 


. 86247 


0.2498 


. 86256 


-.00009 


0.30 


. 83961 


0.2998 


. 83972 


-.00013 


0.35 


. 81796 


0.3497 


. 81810 


-.00017 


0.40 


. 79742 


0.3996 


.79760 


-.00023 


0.477 


. 76766 


0.4765 


.76790 


-.00030 


0.50 


. 75929 


0.4997 


. 75955 


-.00034 



2 

Slope of Gp; versus q using Gross' correction is 
" 0.0218 ± 0.0022. 

Slope of Gp; versus q2 using Friar's correction is 
" 0.0213 ± 0.0022. 
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TABLE D-3 





COMPARISON OF 


Cg(q^) and 


CE(q"') 


FOR LF5 


2 

q 




c^(q') 




2 

q 


CE(q2') 


2 ' 

Cp(q ) 


0.05 




.96881 




0.0500 


. 96881 


-.00000 


0.10 




.93966 




0.100 


.93968 


-.00002 


0.20 




. 88665 




0.1999 


.88670 


-.00006 


0.25 




. 86244 




0.2498 


. 86252 


-.00010 


0.30 




. 83956 




0.2998 


. 83968 


-.00014 


0.35 




. 81790 




0.3497 


. 81805 


-.00018 


0.40 




.79735 




0.3996 


. 79754 


-.00023 


0.477 




.76759 




0.4765 


. 76782 


-.00031 


0.50 




. 75921 




0.4997 


.75947 


-.00034 


Slope of 




versus 

O.i 


q2 

0220 


using Gross 
± 0.0022. 


' correction is 


Slope of 


Gt, 

-t-n 


versus 

0.1 


q^ 

0214 


using Friar 
± 0.0022. 


's correction is 










TABLE D- 


4 






COMPARISON OF 


CE(q^) AND 


CE(q"’) 


FOR LF15 


2 

q 




CE(q') 




2 

q 


CE(q^') 


9 ' 

, > 

^ 2 
CE(q') 


0.05 




.96900 




0.0500 


. 96900 


-.00000 


0.10 




.94002 




0.1000 


.94003 


-.00002 


0.20 




.88731 




0.1999 


. 88737 


-.00006 


0.25 




. 86324 




0.2498 


.86332 


-.00010 


0.30 




. 84049 




0.2998 


. 84060 


-.00013 


0.35 




. 81895 




0.3497 


. 81909 


-.00018 


0.40 




.79851 




0.3996 


. 79869 


-.00023 


0.477 




.76889 




0.4765 


.76913 


-.00030 


0.50 




.76056 




0.4997 


.76082 


-.00034 


Slope of 
Slope of 


Ge 

Ge 


versus 

0. 

versus 

0. 


q^ 

0185 

q2 

0180 


using Gross 
± 0.0022. 
using Friar 
+ 0.0022. 


' correction is 
's correction is 
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APPENDIX E 



COMPUTER PROGRAMS 

Six FORTRAN programs were used to perform calculations 
needed for this work. They are TOPRADCR, TOPROTON, TOPDEUTR, 
BERHFORM, TRADFORM and BERRELCR. The first three were 
written by Topping and can be found in Reference [5] . 



BERHFORM 


and TRADFORM calculate the experimental cross 


sections 


for each of the four channels considered for 


hydrogen 


and deuterium respectively. BERHFORM is repro- 


duced on 


the following pages. BERRELCR was used to 



compare the relativistic corrections of Gross and Friar. 
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